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Abstract. The electrostatic decay enables energy transfer from a finite amplitude Langmuir to a backscattered daughter 
Langmuir wave and ion acoustic density fluctuations. This mechanism is thought to be a first step for the generation of type 
III solar radio emissions at twice the plasma frequency. The electrostatic decay is here investigated through Vlasov-Poisson 
simulations by considering Langmuir localized wave packets in the case T e = T p . Simulation results are found to be in good 
agreement with recently reported observations from the STEREO mission of the electrostatic decay of beam-driven Langmuir 
waves during a type III burst. 
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INTRODUCTION 

Solar Type III radio emissions are one of the most promi- 
nent features of the meter-decameter ranges of frequency. 
The emissions show a pronounced drift with time to- 
wards lower frequencies. Since the early work of Wild 
[1] and Ginzburg and Zheleznyakov [2], the generally 
accepted model for such emission is as summarized be- 
low. During a flare, high energy electrons (1-100 keV) 
are expelled from the solar corona and travel along the 
interplanetary magnetic field lines. The supra-thermal 
electrons produce a bump on the local electron distri- 
bution function generating Langmuir waves via the so- 
called "bump-on-tail instability". Then, nonlinear wave 
couplings generate electromagnetic waves at /~ (the lo- 
cal electron plasma frequency) or 2f~ . The plasma fre- 
quency decreases with the heliocentric distance due to 
the decrease of the electron density: this is the origin 
of the time frequency drift characteristic of the Type III 
emissions. 

Type III electromagnetic emissions are thought to be pro- 
duced via two different nonlinear wave-wave couplings. 
Through electromagnetic coupling, a mother Langmuir 
wave L decays into a low frequency LF waves and a 
transverse electromagnetic wave T f - at the local plasma 

J p 

frequency, observed as Type III fundamental emission: 

L -> T f - +LF 

jp 

Through electrostatic coupling, known as Langmuir 
Electrostatic Decay (hereafter LED), the mother Lang- 
muir wave L decays into a low frequency ion acoustic 
wave S and a daughter Langmuir wave L', which can fur- 
ther non linearly couple with the pump wave to generate 



a transverse electromagnetic wave T 2 j~ at twice the local 
plasma frequency, observed as Type III harmonic emis- 
sion: 

L -t L'+S then L'+L -> T 7f - 

We hereafter concentrate on the LED. 
Henri et al. [3] recently reported direct observations 
of Langmuir waves decaying into secondary Langmuir 
waves and acoustic waves during a Type III solar event, 
from STEREO/WAVES data. They found that: 

■ the Doppler-shifted frequencies of the three ob- 
served waves satisfy the resonant relations of 
momentum and energy conservation expected for 
three-wave coupling 

co L = at L , +(o s k L = k L > +k s (1) 

■ a bicoherence analysis confirms the phase coher- 
ence of the three waves; 

■ the coupling regions are spatially localized with size 
of about 2000 Debye lengths. 

In this former work, the LED threshold and the growth 
rate of IAW density fluctuations generated by LED were 
both evaluated from analytical solutions involving a 
purely monochromatic three-wave coupling [4]. How- 
ever, observations show that: (i) the Langmuir waves are 
isolated wave packets with a packet width of the order of 
a few wavelengths; (ii) proton and electron temperatures 
are known to be close in the solar wind so that IAW asso- 
ciated with the LED should be strongly Landau-damped, 
this would limit the development of the IAW and thus the 
LED. 



The goal of this paper is to study the dynamic of 
the LED through ID- IV Vlasov-Poisson simulations by 
considering an initial localized Langmuir wave packet, 
for equal electron and ion temperature. For a better com- 
parison with observed waveform during Type III bursts, 
the simulation results are also presented as they would 
appear if observed by spacecraft instruments. 



DESCRIPTION OF THE MODEL 

LED is here investigated through 1D-1V Vlasov-Poisson 
simulations in the electrostatic approximation. The de- 
scription of the code and the numerical scheme can be 
found in Mangeney et al. [5]. 

The Vlasov-Poisson system is solved for the 1D-1V 
electron and proton distribution function, f e (x,v,t) and 
f p (x,u,t). The equations are normalized by using the 
following characteristic electron quantities: the charge 
e, the electron mass m e , the electron density n e , the 
plasma (angular) frequency co pe = \J '\%n e e 2 /m e , the De- 
bye length Xd = \jT e /4nn e e 2 , the electron thermal ve- 
locity v t h je = XoCOpe = y/T e /m e and an electric field 
E = m e v t h e (0 P e/e. Then, the dimensionless equations for 
each species read: 
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where v (resp. u) is the electron (resp. ion) velocity nor- 
malized to the electron thermal velocity. Furthermore 
jU = m p /m e = 1836 is the proton-to-electron mass ratio 
and <j) and E are the self-consistent electric potential and 
electric field generated by the plasma charge density fluc- 
tuations according to Poisson equation (Eq. 4). Finally, 
E ext is an "external" driver added to the Vlasov equa- 
tion that can be switched on or off during the runs. The 
electron (resp. ion) distribution function is discretized in 
space for < x < L x , with L x = 10000 Xd the total box 
length, for a resolution of dx — Xd- The electron velocity 
grid ranges over —5 < v/v t ^ e < +5, with a resolution of 
dv = 0.04 v t ^ e . (resp. —5 < u/u t h,i < +5, with a resolu- 
tion of du — 0.04 u^ j for the proton velocity grid, where 
u th,i = \jT p /m p is the proton thermal velocity). The tem- 
peratures are chosen to be equal T p = T e . 

Periodic boundary conditions are used in the spatial 
direction. The electron and proton distributions functions 
are initially Maxwellian with respect of velocity, with a 
random noise in density. 



Finally, the initial Langmuir wave packet with the 
desired wavelength Xl and electric field amplitude El 
is generated by an "external" electric field E ext . The 
pump Eext oscillates at the expected Langmuir frequency 
(Ol = (O pe + 3/2{2%Xd / Xl) 2 , with a phase equal to a>it — 
2%x/Xl, and is spatially localized with a gaussian-shaped 
envelop. A Langmuir wave packet propagating only in 
the forward direction is thus excited. The external elec- 
tric field E ext is switched off when the generated Lang- 
muir wave reaches the desired amplitude El, typically 
at t <~ 300 (Op e l , much shorter than the LED timescale. 
The initial Langmuir wave packet then evolves self- 
consistently. Further details on the forcing can be found 
in [9]. 



LANGMUIR ELECTROSTATIC DECAY 
OF A LOCALIZED WAVEPACKET 

In order to compare the simulations results with observa- 
tions of electric waveforms during Type III events, Lang- 
muir wavelength and amplitude are set as indicated by 
solar wind observations. We choose a Langmuir wave 
packet with wavelengths centered on Xl = 200 Xd, a 
packet width A = 2000 Xd, and a maximum initial elec- 
tric field E L = 6 10~ 2 . Since the electric field associated 
with IAW is known to be low, in the following, ion den- 
sity fluctuations will be used as a tracer for IAW. 

LED is observed in Vlasov simulations for a level of 
Langmuir electric amplitude typical of those observed in 
the solar wind during type III bursts. The evolution of 
the Langmuir wave packet is shown in Fig. 1. The elec- 
trostatic decay starts at t ~ 10 4 CO" 1 . The mother and 
daughter Langmuir wave packets can be follow in the 
top left panel that shows the space-time evolution of the 
electric density energy E(x,t) 2 /2. The Langmuir mother 
wave packet propagates towards the right at its group ve- 
locity (black plain line) and emits daughter Langmuir 
wave packets propagating backward at their own group 
velocity (white dashed line). The bottom left panel shows 
the temporal evolution of ion density fluctuations dur- 
ing the decay of the Langmuir wave. The fluctuations 
along the dashed line is a small amplitude artifact of 
the initial forcing, it does not interact with the elec- 
trostatic coupling. Instead, ion density fluctuations gen- 
erated from the Langmuir mother wave packet (along 
the black solid line) propagate forward at the ion sound 
speed (slope of the dashed line) with the expected wave- 
length for the IAW decay product X 1A ~ 150 X D . They 
are the decay product. Note that the IA fluctuations are 
generated locally where the two Langmuir wave pack- 
ets beats. Indeed, IAW density fluctuations are heavily 
Landau damped, since T e = T p , as soon as the waves es- 
cape the area where LED occurs. Therefore, the LED- 
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FIGURE 1. Evolution of a localized Langmuir wave packet with initial normalized amplitude Ei = 6 10~ 2 . Top panel: electric 
density energy E(x,t) 2 /2. Bottom panel: proton density fluctuations. The expected group velocity of the mother (resp. daughter) 
Langmuir wave is shown by the black plain line (resp. white dashed line). The ion sound speed is shown by the black dashed line. 



produced IAW is expected to be observed in the solar 
wind only locally where the mother and daughter Lang- 
muir waves interact. 

In order to compare the simulation results with in- 
situ waveform observations, we mimic the conditions 
of observation onboard a spacecraft that would record, 
on monopole antennas, a decaying Langmuir wave. We 
thus hereafter introduce in the presentation of simula- 
tion results (i) a spacecraft floating potential effect, (ii) 
a Doppler-shift effect. 

First, monopole antenna channel are known to be sen- 
sitive to both the electric field and the proton density 
fluctuations, through the fluctuation of the spacecraft po- 
tential [6, 7, 8]. This means that the observed signal 
is a combination of both electric field and proton den- 
sity fluctuations. To reproduce this effect in the simula- 
tion, we consider an "equivalent signal" s(t) — E(x,t) + 
a n p (x,t), with a the equivalent in the simulation of a 
calibration parameter that gives in the observations the 
ratio of the spacecraft floating potential with respect to 
density fluctuation. This parameter is set to a = 50 as an 
order of magnitude following the prescription of Kellogg 
etal. [8]. 

Second, observed waveforms in the spacecraft frame are 
Doppler-shifted since the plasma is moving at the solar 
wind speed. A Doppler-shift effect is introduced in the 
simulation by considering the "equivalent signal" s(t) at 
position x(t) = xq + V s t where V s is a constant velocity. 
We choose here V s — Vthe as a first approach. 
This "equivalent signal" s(x(t),t), obtained from the sim- 
ulation, represents the signal that would be recorded by 
a spacecraft crossing such a decaying Langmuir wave 
packet. 

For an easy confrontation with observed waveforms, 
simulation results are analyzed and presented in Fig. 2 
the same way observations have been analyzed and pre- 
sented by Henri et al. [3]. Both a Doppler-shift effect 



and the sensitivity to both the electric field and the den- 
sity fluctuations are taken into account. The top left panel 
shows the electric field (grey) and proton density (black 
line) waveforms from the simulation. Both are plotted in 
normalized units, as described in the simulation model. 
The left bottom panel shows the wavelet transform of 
s(x(t),t). The simulations show that (i) beat-like mod- 
ulation at the plasma frequency would be observed, (ii) 
the IAW signal would be seen at lower frequency cen- 
tered where the maximum of the beat-like Langmuir sig- 
nal is observed. This is in full agreement with the obser- 
vations. The Fourier spectrum of s(x(t),t) is shown in the 
right top panel, with a zoom at the plasma frequency in 
the right bottom panel. The frequency of IAW density 
fluctuations as well as the separation between the fre- 
quency peaks of the two Langmuir waves are essentially 
a Doppler-shift effect. The full line represents the total 
Fourier spectrum. Dashed and dotted lines distinguish 
the contribution of the ion density fluctuations and elec- 
tric field respectively. The simulation thus proves that the 
high frequency part of the observed spectrum is domi- 
nated by the response of the antenna to the fluctuations 
of the electric field, whereas the low frequency part is 
dominated by the response of the spacecraft potential to 
density fluctuations. 

The waveform and spectrum obtained from Vlasov 
simulations agrees qualitatively well with observations 
reported in Henri et al. [3]. However, in order to check 
whether the level of density fluctuations is quantitatively 
consistent with previsions from Vlasov simulations, a 
full calibration of the floating potential of the STEREO 
spacecraft is needed in order to get a better evaluation 
of a that fully takes into account the geometry and 
dimension of the spacecraft. 




FIGURE 2. Waveform, spectrum and wavelet transform from simulation results, as they would appear when observed by 
spacecraft instruments when crossing a Langmuir wave under LED. Top left panel: waveform in a moving frame: electric field 
(grey) and proton density (black). Bottom left panel: Morlet wavelet transform of signal s(x(t),t). Top right panel: corresponding 
Fourier spectrum (full line) with respective contribution of the density fluctuations (dashed line) and electric field (dotted line). 
Bottom right panel: zoom on the double peak at the plasma frequency. 



CONCLUSION 

We have reported ID- IV Vlasov-Poisson simulations of 
the Langmuir Electrostatic Decay. The simulations have 
been done in typical solar wind conditions: (i) equal elec- 
tron and proton temperature, in order to take into ac- 
count the strong Landau damping of daughter ion acous- 
tic waves, and (ii) localized wave packets, in order to 
consider the limited interaction time between the mother 
and daughter waves, each propagating at its own group 
velocity. 

Vlasov simulations in typical solar wind conditions (i) 
shows that the observed level of Langmuir electric field 
during type III burst is high enough for Langmuir Elec- 
trostatic Decay to start, (ii) reproduce qualitatively the 
beat-like Langmuir waveforms as well as the spectrum 
observed during type III bursts when obtained from 
monopole antennas. 

In order to fully confirm the interpretation of the ob- 
served Langmuir waveforms in term of the LED of 
type III beam-driven Langmuir waves, an effective 
threshold will be computed from vlasov simulation in 
typical solar wind conditions [9]. 

We have study here a first step for the generation of elec- 
tromagnetic waves like for example Type III emission. 
Further Vlasov-Maxwell simulations must be performed 
to study the next step of the process and check the energy 
transfer from electrostatic to electromagnetic emissions. 
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